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Recent experimental results from the STAR collaboration suggest event-by-event charge sepa- 
ration fluctuations perpendicular to the event plane in non-central heavy-ion collisions. Here we 
present the correlator previously used split into its two component parts to reveal correlations par- 
allel and perpendicular to the event plane. The results arc from a high statistics 200 GcV Au+Au 
collisions data set collected by the STAR experiment at RHIC. We explicitly count units of charge 
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separation from which we find clear evidence for more charge separation fluctuations perpendicular 
than parallel to the event plane. We also employ a modified correlator to study the possible P-even 
background in same and opposite charge correlations. 

PACS numbers: ll.30.Er, 11.30.Qc, 25.75.Ld, 25.75.Nq 



I. INTRODUCTION 

Parity violation represents a preference of handedness 
in nature. It may be violated globally or locally. In 
the global sense, the weak interactions of the standard 
model are parity odd while the strong interactions 
are parity even at vanishing temperature and isospin den- 
sity However, it has been found possible for parity 
to be violated locally in microscopic domains in QCD at 
finite temperature [31 H]. Parity-odd (P-odd) domains 
in QCD are the consequence of topologically non-trivial 
configurations of gauge fields. A particular domain may 
be characterized by its topological charge. States with 
positive and negative topological charge violate parity 
in the opposite sense. Only states with zero topological 
charge conserve parity. The global conservation of par- 
ity in QCD occurs since positive and negative topological 
charge states are equally probable in nature. 

The hot, dense, and deconfined QCD matter produced 
at RHIC is a natural place to study such P-odd domains. 
Recently a hypothesis has been made stating that these 
P-odd domains might be observable in heavy-ion colli- 
sions. The so called Chiral Magnetic Effect (CME) states 
that P-odd domains can interact with the very large 
magnetic fields in non-central collisions yielding charge- 
separation parallel to the system's orbital angular mo- 
mentum [SMj. This can be viewed as the creation of 
an electric dipole moment vector perpendicular to the 
reaction plane (the plane which contains the impact pa- 
rameter and the beam momenta). 

For a given sign of topological charge and magnetic 
field, the sign of the electric dipole moment produced by 
the CME is also fixed (parity violation). However, pos- 
itive and negative topological charges are equally likely 
and cannot be distinguished on an event-by-event basis. 
One therefore expects the CME to instead manifest it- 
self in an experiment as charge-separation fluctuations 
perpendicular to the reaction plane. 

The reader is referred to the previous STAR publica- 
tions on this subject for more descriptions on the basics of 
this analysis [9l [10] . This article complements the previ- 
ous publications in two principle ways. First, we present 
the correlator, (cos((/)q -|- (/)/3 — 24'rp)), split into its in- 
plane and out-plane components (See Eq. [T]). Second, we 
compare the correlator previously used to a modified cor- 
relator. The comparison enables a better understanding 
of the suppression of opposite w.r.t. same charge corre- 
lations measured with {cos{(j)a + 0/3 — 2^1*^?)). 

This article is divided into six sections. In Sec. II we 
describe the STAR experimental setup and data taking 
conditions used in this analysis. In Sec. HI we describe 
the methodology of the analysis including the definitions 



of correlations measured. In Sec. IV we discuss the sys- 
tematic uncertainties which mainly arise due to event 
plane resolution uncertainties in the modified correlator. 
In Sec. V we present our results. Finally in VI we sum- 
marize our results. 



II. EXPERIMENTAL SETUP AND DATA 
TAKING 

In this analysis, 57 million minimum bias events taken 
by the STAR detector [IT] at RHIC during the 2007 
Au-I-Au run at -^/snn — 200 GeV are used. A hadronic 
minimum bias trigger was formed by requiring a spec- 
tator neutron signal above the threshold value in both 
zero-degree calorimeters (ZDC). Two ZDC shower max- 
imum detectors (ZDC-SMD) measure the spectator neu- 
tron spatial distributions. The ZDC-SMDs are located 
in the beam rapidity regions [12]. Charged particles 
were tracked primarily with the STAR Time Projec- 
tion Chamber (TPC). Tracks are retained if their trans- 
verse momentum and pseudorapidity are in the range 
0.15 < pt < 2 GeV/c and \'q\ < 1.0, respectively. Event 
and track cuts are chosen to be the same as in the previ- 
ous STAR publications on this subject [51 [TO]. Centrality 
in this data set is determined from the global tracking of 
charged particles satisfying specific track quality cuts in 
the pseudorapidity region |?7| < 0.5 and with the dis- 
tance of closest approach (DCA) to the primary-vertex 
less than 3 cm [13) . 

III. METHOD OF ANALYSIS 

The correlation function used in our previous publica- 
tion to search for the CME is given by |14) : 

(cOs((/)„ -t- 0/3 - 2*Rp)) 

= (cos(A0q) cos(A0^) — sin(A0a) sin(A0^)) 

= [(wi.aWi./s) + Sin] - [(ai.Qfli,/?) + -Bout]- (1) 

The averaging is done over all particles in an event and 
over all events. and 0/3 are the azimuthal angles of 
particles a and /?, respectively, ^'rp represents the az- 
imuthal angle of the reaction plane and A0 = (0 — 5'rp). 
-Bin and -Bqut represent P-even background processes 
which may or may not cancel, vi and ai are the first 
harmonic coefficients in the Fourier decomposition of the 
azimuthal distribution of particles of a given transverse 
momentum and rapidity: 




cx l-f 2vi^a cos(A0) + 2v2,a cos(2A0) + ... 
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2ai.Q sin(A0) + 2a2.a sin(2A(/)) 



(2) 



We refer to Eq. [T] as the three-point correlator. Since 
the reaction plane is not directly measurable we estimate 
it using event planes. The event planes are calculated 
from the particle distributions themselves, 



^„ = - tan 
n 



Wi sin(n(/>i) 



Y,Wi cos(n0i) 



(3) 



where n is the harmonic and Wi is a weight for each parti- 
cle i in the sum [15] . The weight is chosen to be the px of 
the particle itself when 0.15 < pt < 2 GeV/c to increase 
the event plane resolution. Above 2 GeV/c, the weight 
is set to 2. A first harmonic event plane is obtained 
from the spectator neutron distributions detected in the 
STAR ZDC-SMD. This type of event plane exploits the 
directed flow of spectator neutrons measured at very for- 
ward rapidity. A second harmonic event plane (^'2) is 
obtained by exploiting the large elliptic flow of charged 
hadrons measured at mid-rapidity in the TPC. 

As the CME causes charge separation fluctuations per- 
pendicular to the reaction plane, it is the sine part of 
the three-point correlator which is sensitive to the CME. 
Note that the cosine part serves to establish a reference or 
baseline to the measurement since both parts are equally 
sensitive to backgrounds unrelated to the reaction plane. 
In this article we present measurements of both parts. 

The three-point correlator weights different azimuthal 
regions of charge separation differently, i.e. oppositely 
charged pairs which are emitted azimuthally at 90° from 
the event plane (maximally out-plane) are weighted more 
heavily than those emitted only a few degrees from the 
event plane (minimally out-plane). We wish to modify 
the three-point correlator such that all azimuthal regions 
of charge separation are weighted identically. The modi- 
fication, in particular, allows us to better understand the 
source of the suppression of opposite charge correlations 
seen previously [9l [10] . 

This may be done by first rewriting Eq. [T] as 

(cos(0„ + (/-/J - 2^'Rp)) = 

{{M^MpSa,Sp)m) - {(M^MpS^SphuT) , (4) 

where AI and S stand for the absolute magnitude (0 < 
M < 1) and sign (±1) of the sine or cosine function, 
respectively. IN represents the cosine part of Eqjl] (in- 
plane) and OUT represents the sine part of Eq. Ill (out- 
plane) . 

To study the dependence of expression [4] on M we com- 
pare the correlations obtained to those of a reduced ver- 
sion: 



((S'qS'^)jP^ - {SaSfs) 



out; 



msc. 



(5) 



We refer to Eq. [5] as a modulated sign correlation (msc) . 
The transition from Eq. [4] to Eq. [5] can be seen with 
the following two reductions: {MS) — {M) (S) and 
{MaMfj) — {M) . {SaSp) may be written as sum of 



terms involving flow coefficients of which only the odd 
harmonics contribute. The common coefficient for all 
contributions is (4/7r)^ with a pre-factor of 1/n, where n 
is the order of the harmonic. For this reason we choose 
{Mmf = (Mout)^ = (7^/4)^- With this choice, the msc 
is also given by the far right hand side of Eq. [TJwhen the 
n = 1 Fourier coefhcients dominate over the other odd 
coefficients. The msc differs from the three-point cor- 
relator in the inclusion of higher harmonics and in the 
removal of " magnitude correlations" (correlations of Mq, 
with M/j), which are of non-flow origin. Since the msc 
is not a pure harmonic, its event plane resolution correc- 
tion is also not generally localized within one harmonic. 
However, we are justified in using the same correction so 
long as Oi ^ On/n and vi ^ Vn/n {n = 3, 5, 7...) or if at 
least a„ fluctuations are similar in magnitude as u„ fluc- 
tuations. We correct both the three-point correlator as 

well as the msc with a second harmonic sub-event plane 

1/2 

resolution, (cos(2(\E'a — ^f,))) ' , where ^'a and vpf, are 
the event plane angles in sub-event a and 6, respectively. 
We discuss the systematic uncertainties associated with 
this correction applied to the msc in Sec. IV. 

For a known reaction plane, (SaSp) is given simply by 
the net number of particle pairing combinations divided 
by the total number of combinations. The net number of 
particle pairing combinations is defined as the difference 
in the number of same side and opposite side combina- 
tions. For the same charge in-plane correlations we have 



{SaSp] 



{N^ - 1) 



iVR (iVR _ 1) _ 2N^N^ 



IN 



1 



(6) 



For opposite charge in-plane correlations we have 



(7) 



where stands for the number of particles detected ei- 
ther on the left (L) or right (R) of the perpendicular 
to the reaction plane in the transverse plane and with a 
positive (-I-) or negative (-) charge. For out-plane correla- 
tions one simply replaces L and R with T and B (top and 
bottom of the reaction plane in the transverse plane). 

To avoid self correlations where a particle is trivially 
correlated with an event plane calculated in the same 
particle pool, we use two equal multiplicity sub-events to 
calculate the msc. The sub-events are statistically inde- 
pendent with random particle assignments. With sub- 
events, the azimuthal locations (T,B,L,R) of the parti- 
cles from one sub-event are calculated with respect to the 
sub-event plane from the other sub-event. 



A. Charge separation counting 

Units of in-plane and out-plane charge separation are 
defined as 

AQiN = - A^L) - (A« - N^) , 
AQouT = {Nj AT) - (atB _ ^B) 
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respectively. They can be understood as the net charge 
on one side of the event plane minus the net charge on 
the opposite side of the event plane. The choice of sign 
for AQ is irrelevant here. 

Equation [s] like Eq. [l] is sensitive to P-even correla- 
tions not related to charge separation. With the aim of 
separating out the simplest effects of charge separation 
fluctuations from other P-even backgrounds we express 
the msc in terms of states of observed charge separation 
AQ. By simplest we mean the contribution to Eq. [5] 
which arises from different AQqut and AQjn probabil- 
ity distributions. We rearrange the msc into two terms: 

msc = Amsc + AiV (9) 
A msc = 

]4 Eaq (^(AQ)> [msciN(AQ) - mscouT(AQ)](10) 
AN = 

^ Eaq (msc(AQ)) [7Vin(AQ) - 7Vout(AQ)] , (11) 

where the sum goes over all observed units of 
charge-separation. N-e stands for the total number of 
events. iViN(A(5) stands for the number of events 
with AQ units of in-plane charge separation. msc(A(5) 
stands for the msc in the bin AQ. The aver- 
ages, {N{AQ)) = (^in(AQ) + 7Vout(AQ))/2 and 
(msc(A(5)) — (msciN(A(5) -|- mscouT(A(5))/2, represent 
an average over in-plane and out-plane parts. 

A given AQ state will be a superposition of many dif- 
ferent configurations or sub-states. The sub-states may 
be described in terms of an underlying neutral pairing 
of particles plus the residual net charge on each side (T, 
B, L or R). The underlying neutral pairs are formed by 
pairing up positively and negatively charged particles on 
a particular side until only a residual net charge remains. 
The residual net charge in a given AQ bin may also be 
arranged in several ways. For example, consider the state 
AQiN = +2. The residual net charge is +2 units. One 
sub-state is the case when the left side has a net charge 
of +2 and the right side has a net charge of (neutral) . 
Another sub-state is formed with a net charge of -2 on 
the right and on the left. The other sub-state occurs 
when the left side has -1-1 and the right side has -1 units 
of net charge. The idea of charge separation counting is 
illustrated in Fig. [l] Both the underlying neutral pair- 
ing and the residual net charge contribute to the overall 
configuration within each AQ bin. 

The right hand side of Eq. [9] is composed of two 
terms. The first term, Amsc, is sensitive to the differ- 
ence between in-plane and out-plane AQ configurations 
(msciN(AQ) — mscouT(AQ)). The second term, AA^, 
is sensitive to the difference between in-plane and out- 
plane AQ probabilities {Nw{AQ) - AfouT(AQ)). The 
factor, Ain(AQ) — Aout(^Q), is of course identical 
for same and opposite charge correlations. Therefore, 
the difference between same and opposite charge AA^ 
correlations is determined exclusively by the prefactor 
(msc(AQ)). If the CME does not significantly alter the 
charge-separation sub-states it will be isolated in the AA^ 
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FIG. 1: (Color online) Example of a charge configuration 
with the underlying neutral pairing enclosed by dotted ovals, 
(a) shows the procedure for counting AQout ~ +5. (b) 
shows the same event but with the procedure for counting 
AQiN = +1. 



term of Eq. [9] In general, both terms could be affected 
by a P-even background and neither is to be regarded as 
an isolation of a P-even background. 

The effects of P-even local charge conservation coupled 
with non-zero V2 has been shown to yield a substantial 
background to the same charge correlations of Eq. [l] [TB] . 
The contribution of this effect to Eq. [5] should be reduced 
due to the treatment of the magnitudes of the cosine and 
sine functions. 

The effects of momentum conservation, in general, will 
be equally present perpendicular and parallel to the reac- 
tion plane. The focusing of particles closer to the reaction 
plane, V2 > 0, which causes these effects to be enhanced 
in the cosine parts relative to the sine parts in Eq. [TJ is 
effectively undone when we now consider Eq. [5j 



B. Acceptance effects 

Anisotropic or imperfect detector acceptance may also 
induce false correlations. The STAR detector has nearly 
complete azimuthal coverage but nevertheless we apply 
a re-centering correction [17 to all the event plane calcu- 
lations. The correction is done in bins of centrality, loca- 
tion of collision parallel to the beam axis, and STAR run 
number which represents a period of time with constant 
detector calibrations. Particles a and /3 in the three-point 
correlator are also re-centered. The effect of this proce- 
dure was only found to be sizable in the most central bins 
where the signal is small. 
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IV. SYSTEMATIC UNCERTAINTIES 

We estimate the systematic uncertainties on our mea- 
surements by comparing results obtained using TPC and 
ZDC-SMD event planes. The difference between these 
two measurements divided by the statistical error forms 
our estimate of the systematic uncertainty. This estimate 
is shown in the shaded bands in Figs.[5j[6j and[7] For the 
three-point correlator, these values characterize non-flow 
uncertainties in the reaction plane reconstruction. For 
the AA^ and Amsc terms the values characterize both 
non-flow uncertainties in the reaction plane reconstruc- 
tion as well as uncertainties in applying the 2"^^ harmonic 
event plane resolution to the msc. 

Other systematic uncertainties were studied exten- 
sively in the previous publications on this subject p,[TU]. 
All were shown to be negligible compared to the uncer- 
tainty in determining the reaction plane. The shaded 
bands in the figures here represent the same uncertainty 
determined by a comparison of measurements with f** 
and 2"'^ harmonic event planes. 

For the simplified case of pure elliptic flow {v2 > 0) -I- 
an added CME signal (|ai| > 0), we have also verified 
through Monte Carlo simulations that the msc with re- 
spect to the sub-event planes corrected by the sub-event 
plane resolution is equivalent to the msc with respect to 
the reaction plane. 



V. RESULTS 

Larger charge separation fluctuations perpendicular to 
rather than parallel to the event plane can be seen by 
comparing distributions of AQout to AQin as shown 
in Fig. [2] for the 40 - 50% centrality bin. Figure [s] 
shows the difference over the mean of the RMS values 
( (RTsorT+RS)/2 ) ^^'^^^ centrahty. The CME will 
cause wider out-plane distributions, however, P-even pro- 
cesses may also cause the same feature (e.g. the decays of 
elliptically flowing resonances). Figures [2]and|3] are not 
corrected for the event plane resolution, however they 
clearly demonstrate larger charge separation fluctuations 
perpendicular rather than parallel to the event plane. 

Figure[4]presents (sin(0c — VPi)) for positive and nega- 
tive charges. Such a measure is sensitive to global parity 
violation of the strong interactions, i.e. a preference of 
charge separation orientation relative to the angular mo- 
mentum orientation of the system. The results of Fig. |4] 
do not show a significant charge dependence. The mean 
values of both positive and negative charges are less than 
5 X 10^* at the 95% confidence level. For the most cen- 
tral and peripheral collisions we observe non-zero values 
for ai. However, the values have the same sign for both 
charge types which is inconsistent with a global violation 
of parity. 

The three-point correlator measured with T'* and 2"^^ 
harmonic event planes is shown in Fig. [5j We find con- 
sistency between correlations obtained with both event 




-40 -30 -20 -10 10 20 30 40 

A Q 

FIG. 2: Sample AQ distributions for the 40-50% centrality 
Au-|-Au collisions at ,/snn= 200 GeV. Not corrected for event 
plane resolution. 




80 70 60 50 40 30 20 10 

% Most Central 

FIG. 3: (Color online) ARMS/ (RMS) versus centrality for 
Au+Au collisions at ^/smi= 200 GeV. Not corrected for event 
plane resolution. Errors are statistical only. 



plane types. As the pseudorapidity gap between the 
ZDC-SMD(*i) and the TPC(particles a and (3) is rather 
large (^ 7 units in rj) , we find "direct" three-particle 
effects (clusters) to be an unlikely source for the sig- 
nal. This is an indication that the signal is likely a gen- 
uine correlation with respect to the reaction plane. Also 

5] are our previous results 
which are consistent with 



shown for comparison in Fig, 
from the 2004 RHIC run P 
the current results within statistical errors. 

The modulated sign correlations are compared with 
the three-point correlator in Fig. |6] It is evident that 
the msc is able to reproduce the same trend as the 
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FIG. 4: (Color online) {sm{tfia — ^Pi)) for positive and nega- 
tive charges versus centrality for Au+Au collisions at ^/sm^^= 
200 GeV. Shaded areas represent the systematic uncertainty 
obtained by comparing correlations from positive and nega- 
tive pseudorapidity. 



three-point correlator although their magnitudes differ 
slightly. It is also clear that the correlation magnitude 
for same charge pairs is larger than for opposite charge 
pairs for both correlators. Before any possible interac- 
tion with the medium, the CME is expected to gener- 
ate equal correlation magnitudes for same and opposite 
charge pairs. It was previously supposed that medium 
suppression of back-to-back phenomena could be respon- 
sible for this magnitude asymmetry [9j [10] . Oppositely 
charged pairs from the CME may not freeze out back- 
to-back, but instead with one of the particles deflected 
closer to the event plane due to multiple scattering within 
the medium. This is most likely to occur for the particle 
traversing the largest path length through the medium. 
However, when we weight all azimuthal regions of charge 
separation equally, as with the msc in Fig. [6j we do not 
recover a magnitude symmetry. 

The two terms of the msc in Eq. [9] are shown in Fig. [7] 
We observe that same and opposite charge correlations 
in the AN term have very similar magnitudes, but oppo- 
site signs for all centrality bins. This feature is expected 
from the construction of the AN term due to the rela- 
tively large and approximately equal positive and nega- 
tive charge multiplicities. A model calculation including 
statistical-fdynamical fluctuations of particle azimuthal 
distributions should be performed in order to rule out 
P-even explanations. The Anise term is typically nega- 
tive and has a similar magnitude for same and opposite 
charge correlations. Thus, the source of the magnitude 
asymmetry between same and opposite charge correla- 
tions about zero as shown in Fig. |6] is isolated in the 
Amsc term (Note that the sum of both terms yields the 
total msc). 
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FIG. 5: (Color online) Three-point correlator, Eq. [l] mea- 
sured with 1"* and 2""^ harmonic event planes versus centrality 
for Au-|-Au collisions at y/awN= 200 GeV. Shown with crosses 
are our previous results from the 2004 RHIC run (Y4) 9, IC^. 
The Y4 run used a second harmonic event plane. Y4 and *I'2 
results are consistent within statistical errors. Shaded areas 
for the 2"^^ harmonic points represent the systematic uncer- 
tainty of the event plane determination. Systematic uncer- 
tainties for the V^^ harmonic points are negligible compared 
to the statistical ones shown. 



We now present the composite parts of the three-point 
correlation, Eq. [T] differentially versus 77 andp-r- Figure 
|8] presents the three-point correlator versus the average 
■q of particles a and /? ((77)) and absolute value of the 
difference (IA77I). Figure [9] shows the same composite 
parts versus (pr) and ApT- 

The subtraction of out-plane from in-plane compos- 
ite parts yields the original three-point correlator while 
the sum yields a two particle correlation, (cos((/)c, — 4>p)). 
The split correlations reveal the underlying P-even back- 
ground affecting both composite parts as each part is 
sensitive to event plane independent correlations. We 
see that in each case the functional shape of in-plane and 
out-plane parts are similar. The magnitudes of in-plane 
and out-plane parts are more different for same charge 
pairs. 

Femtoscopic correlations at low relative momentum are 
visible in Figs. [8^-|9]d. The sharp increase of the corre- 
lation strengths for the lowest bins in Figs. [Sj^, |9^, and 
[9]d are due to the combination of quantum interference 
("HBT") in the same charge channel and the final-state- 
interactions (Coulomb dominated) in both channels. Low 
relative momentum in the transverse plane is clearly best 
visible for in Fig. [9]d low values of |Ap7-|. The same phe- 
nomena are also visible for low values of {pt) since these 
values best isolate low values of |Apt|- Low relative mo- 
mentum along the beam axis is clearly visible in Fig. [8]d 
for low values of IA77I. The same phenomena are only 
visible for the larger values of (77) in Fig. [sk since 77 is 
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FIG. 6: (Color online) Modulated sign correlations (msc) 
compared to the three-point correlator versus centrality for 
Au+Au collisions at ^/smi— 200 GeV. Shown with triangles 
is the msc, Eq. Jsl Circles and squares show the three-point 
correlator, Eq. ^ Shaded areas represent the systematic un- 
certainty of the event plane determination. 
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FIG. 7: (Color online) The msc split into 2 composite parts 
versus centrality for Au-|-Au collisions at -^Snn^ 200 GeV. 
Shown with triangles is the Amsc term. Circles and squares 
show the AN term. Shaded areas represent the systematic 
uncertainty of the event plane determination. 
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FIG. 8: (Color online) Three-point correlations split 
up into out-plane (sin(A(;/)a) sin(A(;i/3)) and in-plane 
(cos(A0c«) cos(A(j!)/3)) composite parts for Au-|-Au colli- 
sions at ^/stiN— 200 GeV. (a) shows the correlations versus 
(rj) — {rja + rifi)/2. (b) shows the correlations versus 
= ^ "UfA- Statistical errors are smaller than the 
symbol size. Systematic errors are given by the shaded bands 
and apply only to the difference of out and in plane parts. 



signed. That is, the lowest values of {rj) contain a sub- 
stantial fraction of pairs with the opposite sign of 77 and 
therefore large relative momentum along the beam axis. 

We also observe that the positive signal for oppo- 
site charge correlations observed in the peripheral bins 
of Fig. 15] (cos(A0q) cos(A(/)^) — sin(A(/)Q,) sin(A(/)^)) is 
largely found in the femtoscopic regions of Figs. [S^-|9]d. 
Femtoscopic correlations are sensitive to the size of the 



emission volume at freeze-out [THl HI]- The difference 
between in-plane and out-plane correlations in the fem- 
toscopic region can be due to a difference in the emis- 
sion volumes probed by in and out-plane parts. Such 
a difference may arrise from an azimuthally anisotropic 
freeze-out distribution coupled with flow. 
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charge separation are weighted equally. The same quali- 
tative signal was found to persist in this scheme as well. 
The signal shown in Fig. [6] is largely determined by the 
sign (±) of the cosine and sine functions in Eq. [T] 

We also explicitly counted units of charge separation 
with which we could better understand the source of the 
opposite charge suppression. A parity conserving back- 
ground is more likely to explain the suppression rather 
than the medium induced back-to-back suppression pre- 
viously supposed [nUTU] . A comparison of the RMS values 
for AQouT and AQin suggests greater charge separation 
fluctuations perpendicular to rather than parallel to the 
event plane. The CME as well as P-even processes such 
as flowing resonance decays may both contribute to this 
feature. 

The differential analysis of the in-plane and out-plane 
parts of the three-point correlator versus rj and pt re- 
veals femtoscopic contributions at low relative momen- 
tum. The positive signal in Fig. [s] for opposite charge 
correlations in peripheral collisions is largely found in the 
low relative momentum regions of Fig. [8] and Fig. [9j This 
can possibly be explained by the final-state-interactions 
(P-even) of different emission volumes probed by in-plane 
and out-plane parts. 

Excluding low relative momentum pairs will signif- 
icantly reduce the positive contributions to opposite 
charge correlations in Fig. [Sj However, the difference be- 
tween same and opposite charge correlations will remain 
largely unchanged and consistent with the expectations 
of the CME. P-even local charge conservation coupled 
to elliptic flow modeled by charge balance functions has 
also been shown to generate same charge three-point cor- 
relations comparable to the observed one [16]. A careful 
calculation of the mentioned P-even backgrounds needs 
to be made before a further assessment of the Chiral Mag- 
netic Effect can be made in heavy-ion collisions. 



FIG. 9: (Color online) Three-point correlations split up into 
out-plane and in-plane composite parts for Au-f Au collisions 
at ysNN= 200 GeV. (a) shows the correlations versus (pr) = 
{pT,ci +Pt,;3)/2. (b) shows the correlations versus |Apt| = 
\PT,a — Pt,^|. Statistical errors are smaller than the symbol 
size. Systematic errors are given by the shaded bands and 
apply only to the difference of out and in plane parts. 



VI. SUMMARY 

Correlations sensitive to charge separation in heavy- 
ion collisions have been presented. Consistency between 
correlations with respect to 1** and 2"'^ harmonic event 
planes demonstrates that the signal is likely to be related 
to the reaction plane. Also presented was a reduced ver- 
sion of the three-point correlation in which all regions of 
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